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Last years many scientists used symplectic integrators [1], which are much faster than usual 
integrators. For asteroids, a step of integration for a symplectic integrator usually was taken in 
a range between 7 and 30 days. Some scientists compared the results obtained with different 
> ■ integration steps, but usually they did not compare these results with those obtained with a 
usual integrator. We have made series of computer runs of the orbital evolution of asteroids 
ff~j . for the 3/1 and 5/2 resonances with Jupiter using both the symplectic integrator RMVS3 [1] 
and the Bulirsh-Stoer method [2] (BULSTO). The results obtained with different integrators 
| and different integration steps d s were compared in order to understand what error we usually 
make if we use a symplectic method for investigations of orbital evolution of resonant main-belt 
^ ■ asteroids. For BULSTO the error e per integration step was taken to be less than 10 -8 or 10~ 9 . 

For RMVS3, we have made integrations with d s equal to 3, 10, and 30 days. In each run we 
considered the Sun, 7 planets (except Mercury and Pluto) and N asteroids moving in the 3/1 
or 5/2 resonances with Jupiter (a =2.5 or a =2.823 AU). Initial eccentricities and inclinations 
were the same in all runs: e o =0.15 and « o =10°. Initial values of the mean anomaly and the 
longitude of the ascending node were different. The considered time interval T$ equaled several 
Myr. 

Using orbital elements obtained with a step equal to 500 yr, we calculated the probabilities 
of collisions of asteroids with the terrestrial planets and obtained (for all time intervals and 
all bodies) the total probability Ps of collisions with a planet and the total time interval Ts 
during which perihelion distance q of asteroids was less than a semimajor axis of the planet. 
The values of P r =10 6 P=10 6 P^/N and T=Ty,/N are presented in the Table together with the 
ratio r of the total time interval when orbits were of Apollo type (at a>l AU, q=a(l — e)<1.017 
AU, e<0.999) to that of Amor type (1.017<g<1.33 AU); r 2 is the same as r but for Apollo 
objects with e<0.9. 
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Table: Values of T (in kyr), P r =10 6 P, r, r2, and r^ c for the terrestrial planets (Venus=V, 
Earth=E, Mars=M) at iV=144, Ts=10 Myr (except for the first line for each resonance, for 
which T s =50 Myr). 
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For the asteroids initially located at the 3/1 resonance with Jupiter, we found that the ratio 
of the number of asteroids ejected into hyperbolic orbits to that collided with the Sun is 
much larger for BULSTO than for RMVS3. Besides the values of r^ c presented in the Table at 
7V=144, for the 3/1 resonance at iV=24 we obtained r hc equal to 4.0, 1.7, 0.33, 0.4, and 0.7 at 
£=10~ 8 , 10~ 9 , d s — 3, 10, and 30 days, respectively. So in some cases a symplectic method can 
give large errors. For the 5/2 resonance, the ratio of the values of r^ c obtained by BULSTO 
and RMVS3 also was not small (>3). The difference in values of T and P r was not considerable 
for RMVS3 at d s =10 days and for BULSTO. For d s =30 days at the 5/2 resonance, 78% of the 
probability of collisions with the Earth were caused by 3 asteroids (64%, by two asteroids) and 
52% of all collisions with the Earth were from Aten orbits. 
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